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I. Introduction 


The development of more efficient gas turbine engines and power plants lor future super- 
sonic transports like the high speed civil transport (1ISCT) depends upon the advancement 
of new high temperature materials with temperature capabilities exceeding those of current 
nickel base superalloys. Ordered intcrmctallic alloys, in particular the 11*2 structured NiAL 
have long been considered a prime candidate to replace superalloys in the combustion and 
turbine sections of aircraft gas turbine engines due to a number of property advantages 
including the development of alloying schemes for enhanced creep strength [l]. 

One of the most basic of these strategies is to combine single crystal processing and al- 
loying wit h reactive elements such as Ti. Ilf, Zr and ia in order to generate extremely creep 
resistant materials [2], This approach has been so effective that single crystal NiAl alloy 
turbine vanes have been successfully engine tested [3], Not only did t he parts survive engine 
testing but they displayed superior performance to superalloy components undergoing the 
same evaluation. However, alloy design, while successful to this point, has been entirely 
empirical in nature. Furthermore, limited understanding oi the basic microstruct tires and 
structure-property relationships has occurred due to this empirical alloy design approach. 
For example. Ti additions on the order of 2. 5-3.0 at. % result in a 200-5000 fold reduction in 
creep rate over binary NiAl single crystals [4], The mechanism behind this largo solid solu- 
tion strengthening phenomenon is hampered by the lack of detailed struct ural information 
concerning the NiAI-Ti system. 

While it is well known that ternary and higher order alloying additions to structural 
intermetallics are essential for the optimization ol physical, chemical, structural, and mo- 
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chanical properties. it is only recently that substantial theoretical work has been directed 
towards a thorough understanding of the atomic processes involved, with the goal of elu- 
cidating the role of such alloying additions in controlling properties. In the simple case of 
1VAI and NiAl, then' are abundant, sometimes conflicting results from a variety of tech- 
niques. providing some guidance to the defect structures in these binary compounds as a 
function of stoichiometry [1.5-7]. From a theoretical standpoint, first- -principles calculations 
have also provided very valuable information regarding this issue [s]. However, a full un- 
derstanding of the processes involved and a complete description of the composition and 
temperature dependence of defects and site distributions are not yet available for NiAl and 
I eAl. This problem is even more pronounced for alloy systems which have not been as thor- 
oughly studied and is nearly intractable for ternary and quaternary additions loan ordered 
AH system. The uncertainties and shortcomings of the techniques - whether theoretical 
or experimental used to characterize the behavior of ternary and higher order additions 
raise questions as to whether the problem of highly alloyed systems is easily amenable to 
analysis. 

While it would be extremely advantageous to perform accurate quantital ive and qual- 
itative theoretical analysis of these multicomponent systems, this is beyond the scope of 
most current available techniques. Theoretical predictions of alloy properties are generally 
obtained by two approaches, first -principles and semiempirical. First -principles methods 
involve solving Schrcidinger s equation for the system at hand, while semiempirical methods 
attempt to develop approaches to model energetics with empirical input to obtain unknown 
parameters. Ideally, first principles calculations are best suited for providing the most ac- 
curate and consistent framework for such studies. However, t he small differences in energy 
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related to basic issues like preferential site substitution impose the need for intensive elec- 
tronic structure calculations. Moreover, the complexity of the multi-component systems 
imposes further limitations, and thus, not enough information can be realistically expected 
from a first-principles approach because t he CPF intensive calculations must be performed 
for many possible, large configurations. This limits the usefulness of first principles ap- 
proaches as economical predictive tools. 

On the other hand, semiempirical approaches partially solve this problem by introduc- 
ing a much needed degree of efficiency, which enables one to look at larger systems, at the 
expense of a hopefully minimal loss of physical accuracy. Most importantly, they are useful 
in providing a global view of the mechanisms or processes at hand. However, their range 
is limited, in most cases, to a few elements or some particular crystallographic structure. 
Methods that exhibit a great deal of accuracy for bulk materials fail to reproduce some of the 
most basic properties of surfaces or grain boundaries. In addition, their field of applicability 
to monatomic systems or very specific binary systems (for which specific parameterizat ions 
or potentials are often developed) is limited due to the lack of transferrability of the pa- 
rameters used. Recently, a new semiempirical technique was developed which avoids the 
current limitations on existing methods. The BFS (Bozzolo- Ferrante-Sniith ) method for 
alloys [9] is particularly designed to deal with complex systems and geometries; it has no 
constraints regarding the atomic species under consideration, their number, or the crystal- 
lographic structure. The BFS method provides a simple framework for largo scale computer 
simulations with the appropriate formal background, without the limitations that exist in 
competing techniques. 

After a brief introduction of the BFS method, this paper concentrates on its application 
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to the study of Ni-Al-Ti alloys and the fundamental issues concerning physical properties. 


resulting defect structures, the solubility of Ti iu NiAl and precipitation of a second phase. 
Moreover, the results include Monte Carlo temperature-dependent large-scale simulations 
providing some insight on the relationship between heat treatment and microstructure. 
Finally, owing to the engineering significance of NiAl-Ti alloys and to develop confidence in 
our modelling abilities, the theoretical results derived from the BFS method are compared 
to the results from a concurrent transmission electron microscopy study. 

II. The BFS Method 

Since its inception a few years ago. the BFS method has been successfully applied 
to a variety of problems [9-19]. starting with the basic analysis of bulk properties of solid 
solutions of fee and bcc binary alloys (heat of formation [10-11]. lattice parameter [12]. etc.) 
to more specific applications like the energetics of bimetallic tip-sample interactions in an 
atomic force microscope [id], and Monte Carlo simulations of the temperature dependence 
of surface segregation profiles in Cu-Ni alloys [Id] . An additional advantage of the BFS 
method is that it allows for deriving simple, approximate expressions which describe the 
t rends in segregation and elucidate the driving mechanisms for these phenomena [IS], More 
recently, the ability of the BFS method to deal with alloy surfaces as well as phase stability 
in more complicated systems has been successfully demonstrated in simulations of surface 
alloys containing as many as four elements [17] and the design of Ni-based quaternary 
alloys [9], providing a strong foundation for the work presented in this paper. Also, as a 
consequence of the ideas underlying the foundation of BFS. simple expressions for predicting 
the composition dependence of bulk alloy properties have been recently derived (the BF rule) 
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[19], based solely on pure component properties, providing an alternative lot lie commonly 
used Yegard’s law [20]. 

In wliat follows we present- an operational review of the method. Due to its novel wav of 
partitioning the energy in different contributions, this presentation should be complemented 
with a review of previous applications [9- IS], in order to familiarize the reader with the main 
concepts discussed below. 

The BI S method is based on the idea that the energy of formation of an arbitrary 
alloy structure is the superposition of individual cont ributions s, of nonequivalent atoms in 
the alloy [9], 


= S* + fl ,( f ■' 


( 1 ) 


so that the total energy of formation is 


a// - yy, (2) 

i 

For each atom, we break up the energy into two parts: a strain energy r; s and a chemical 
contribution s chrm . Ihe chemical energy consists of two terms: one {rar( psUon] which 
takes into account the actual chemical environment seen by atom i and a reference term ^ °. 
introduced in order to ensure that no structural information is introduced in t he calculation 
of the chemical effects: z'- hf m = s ( ; - s ( ; 0 , The strain and chemical energy are then linked 
by a cou])Iing factor <j\ 

ti = sf" 1 " +g;s r i hfm 0 ) 

where / denotes t he atomic species of a given atom (c r ° is a reference energy to he defined 
later). 

The strain energy. accounts for the actual geometrical distribution of the atoms 
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surrounding atom i. computed as if all its neighbors were of the sanu atomic species as 
atom /. z* is then evaluated with any suitable technique such as Equivalent Crystal Theory 
(i:<T) [2i]. 

The coupling term, is related to the strain energy in the sense that it contains 
information on the structural defect included in In order to establish this connection, 
based on the assumption that the universal binding energy relationship of Hose et al. [22] 
contains all the relevant information concerning a given single-component system, we can 
write sf as 

( 1 ) 


where 

= !-(! + a m )< ~ (l \ 




and where r/; s ~, given by 




((>) 


is a scaled lattice parameter related to r/f , a quantity that contains the structural infor- 
mation of the defect crystal, a' ./, and are the equilibrium lattice parameter, scaling 
length, and cohesive energy, respectively, of a pure crystal of species i. and q = for 

bcc metals. Once 5; s is evaluated. of* can be easily obtained from Eq.(4) from which the 
coupling term <j t becomes 

9i = (7) 


As in previous efforts [9-19], we chose EOT [21] to perform the strain energy calcu- 
lations. the choice being guided by the simplicity and reliability of this technique. Using 
E(T for computing sf introduces the added advantage that u; s (and thus a ** ) is directly 
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obtained by solving the E(T equation for the defect crystal, as shown below. Within the 
framework of E(’T. a* is interpreted as the lattice parameter of an ideal, perfect crystal 
(i.e.. the equivalent crystal) where the energy per atom is the same as the energy ol atom 
i in the actual, defect crystal. 

In general, the E(T equation for computing the strain energy reads 

A'/fjV -0 h ' ] + = J2 r r j (- {l ' + S{r -> ))r -> (X) 

.1 

where the quant ities p, o. A and t he screening function S arc' defined in Ref. 21 . I lie sum 
on the right hand side of Eq. (X) runs over all neighbors of atom i at a distance r r Lq. 
(S) is then solved for the lattice parameter of the equivalent crystal e; s . R ] and R 2 are 
the corresponding nearest -neighbor and next-nearest-neighbor distances in the equivalent 
crystal. The strain energy is then computed with Eq.(d). For the particular case where 
all the neighboring atoms are located at lattice sites. r t — /q and *S(/q) — 0 for nearest - 
neighbors; r , — r 2 and S( r 2 ) = 1/A for next-nearest-neighbors: and if n is the actual number 
of nearest -neighbors and m is t he corresponding number of next -nearest -neighbors, then Eq. 
(X) is simply 

NR^r~ o}u + + wv! 2 f- {a+ ^ )r \ ( 0 ) 

Rigorously, the computation of the strain energy includes four terms (see Ref. 21). 
In this work, we neglect the three- and four-body terms dealing with the bond angle and 
face-diagonal anisotropies and retain only the two- body term that accounts for bond-length 
anisotropies [21], which we expect to he relevant for atoms in the top (surface) layers. 
The higher order terms would be proportional to the small local fluctuations of the atomic 
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positions around the equilibrium lattice sites. We expert that the leading term of Eq. (4). 


will adequately account for these small distortions. 

The chemical contribution .f; is oblained by an ECT-like calculation. As opposed 
to the strain energy term, the surrounding atoms retain their chemical identity . but are 
forced to be in equilibrium lattice sites. If A’,-* (Mik) denotes the number of nearest ( next )- 
neighbors of species A* of the atom of species i in question, then the f(T equation [JK2I] to 
be solved for the equivalent lattice parameter ei ( ; is 


A /n 




m ii i : 




E 


V, 


+ E^* r 2' 


, — ( li tk + T ~ ) r 2 


(JO) 


where A (M) is the number of nearest (next, (-neighbors in t lie equivalent crystal of species 
i and R\(R 2 ) is the nearest (next (-neighbor distance in the equivalent crystal of lattice 
parameter a ( ; . r \ and r 2 . are the equilibrium nearest- and next-nearest-neighbor distances 
in an equilibrium crystal of species /. respectively. The chemical interaction between atoms 
i and A* is represented by the parameters 077., given by 


<*ik — n i + ^ki ( 1 I ) 

which, extending the underlying concepts in E(’T, parameterizes the tail of the wave func- 
tion in the overlap region between atoms i and A\ The pure element ECT parameter 07 
is then ‘pert urbed' by the BFS parameter A*.,. Therefore, the BFS method introduces 
only two new parameters in addition to the ECT parameters needed for the individual con- 
stituents: A ,7. and A/., for every pair /,/.■ of alloy components. The chemical energy is then 
computed with 
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and 


c-f° = -,o/-r/'*(«r°‘) (1») 

whore ( ~ /t , ) = +1 if > 0 and = -I otherwise, and a] * = q{a l ; - (i‘,)/l;- 

The scaled lattice* parameter r/; " is obtained from lapf 10) and e/-°* is computed by solving 
Kq.( 10) again, but with a;*. = n,. 

The BFS parameters A ^ and A^a in Eq. (11) are obtained from results of first- 
principles. all-elect roll, density-functional based calculations of the* elemental constituents 
and ordered binary compounds of these elements. The particular implementation used in 
this work is the Li near- Muffin- Tin Orbitals (LMTO) method [2d] in the Atomic Sphere 
Approximation (ASA). This scheme was used to calculate the equilibrium properties of the 
elemental solids in the same crystal symmetry as that of the compound to be studied. This 
set of parameters is accurately described by the Local Density Approximation [21]. Thus, 
for this case, we have calculated the properties of bcc-Al. bcc-Ni and bcc-Ti, as well as B2 
ordered NiAL NiTi and TiAl (while the B2 NiAl and NiTi phase do exist in nature. TiAl 
forms a Ll 0 fee- based structure). Calculations wore made for different values of the lattice 
constant, and total energies were then fitted to the universal equation of state of Rose et 
al. [22]. The LMTO method uses a minimal basis set: in this work, we have used only .s, 
p and (I orbitals. All calculations were done with equivalent sampling of the Brillouin zone 
using, for the bcc lattice, 120 A*- points in the irreducible wedge. Apart from the parameters 
describing the equation of state of the element (lattice parameter, cohesive energy and bulk 
modulus), the parameterization of the BFS approach requires the formation energy of a 
single vacancy (E rfI[ .) in order to fix the EOT parameter a. We have also calculated E r „ t . 
with the LMTO method using a supercell approach. Studies of the convergence of this 
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property as a function of the supercell size showed that, for the required precision in the 
calculation ( ~ 0.1 oY/atom) and within practical computational limits, a supercell of S 
atoms is necessary. As BFS is parameterized without considering the relaxation caused by 
the formation of a vacancy, no relaxation is allowed in the LMTO calculations. 

As mentioned above, the consistent parameterization of the BFS method requires also 
the calculation of the formation energy and lattice constant of an ordered binary compound. 
VYe have chosen the bcc based B2 structure. This calculation for the compound is equivalent 
in the basis set and sampling of the Brillouin zone used for the pure elements. 

The EC’T and BFS parameters used in this work for Ni, A1 and Ti are listed in Fables 1 
and 2. respectively. For consistency, these parameters have all been determined purely from 
LMTO calculations and did not involve any experimental input. Once these parameters are 
computed, they remain the same for any calculation involving NiAFTi, requiring no further 
adjustment or replacement . The parameters enter the calculation of the energy of formation 
of the different configurations to be studied through Eqs. ( 4 )-( 10). The numerical values 
of the BFS parameters A ab and A ba can be easily obtained, since one of the advantages 
of the BFS method is that it allows for a simple analytical procedure for the determination 
of such parameters, as described in detail in Appendix A. 

We should emphasize that in the context of BFS. the strain and chemical energy contri- 
butions differ substantially in meaning from the ones these terms have in other approaches. 
The BI S strain energy is related to the usual strain only in that the atomic locations are 
those found in the actual alloy: the BFS strain energy of a given atom is then the actual 
st rain that it would have in a monatomic crystal of the same species of the reference atom. 
Likewise, the BFS chemical contribution is related to the usual chemical energy in t hat the 
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LMTO results 

F(T parameters 

Atom 

Lattice 

( ‘ohesive 

Hulk 

Vacancy 

P 

n 

A 

i 


Pa ram a ter 

Hnergy 

Modulus 

Hnergy 


(A" 1 ) 

(A) 

(A) 


(A) 

(oV) 

(CPa) 

(e\) 





Ni 

2.752 

5.S09 

219.2 

:}.() 

0 

3.0070 

0.703 

0.2710 

Al 

3.102 

3.9-12 

77.5 

l.x 

1 

1 .X750 

l .03X 

0.3005 

Ti 

3.213 

0.270 

121.0 

2.0 

0 

2.0N05 

1.0 FS 

0.372S 


Table I : LMTO results for the lattice parameter, cohesive energy, bulk modulus and vacancy 
formation energy for the bee phases ol Ni. Al and I i. 1 lie last lour columns display the 
ensuing K(T parameters determined from the LMTO results. 


B.I 

.S. Parameters 

All 

A,,h 

A«, 

Ni- Al 

-0.05813 

0.0*22 

Ni-Ti 

-0.005*7 

0.1010 

Al-Ti 

-0.00300 

0.2283 


table 2: HI S parameters A \h and A#i for Ni-AL Ni- Ti and Ti-AI determined by fitting 
the lat tice parameter and heat of formation of the corresponding B2 compounds via LMTO 
calculations. 

actual chemical composition of the alloy is taken into account, but with the neighboring 
atoms located in ideal atomic sites: the BFS chemical energy of a given atom is then the 
actual chemical energy in an ordered environment with the lattice spacing characteristic of 
the equilibrium lattice of the reference atom. Wo refer the reader to previous applications 
of BFS for more insight in this issue [9,19]. 

III. Ni-Al-Ti alloys 

From a practical standpoint two-phase alloys based on a N i A 1 matrix reinforced by 
Heusler (Ni^AlTi) precipitates have been receiving increased experimental attention due to 
their potential as high temperature structural alloys [2-1.25-2*]. The Heusler phase has a 


NASA TM-1 13121 


11 





cubic l.2j structure, in which 1 he unil cell comprises eight bcc unit cells with Al and Ti 
atoms occupying two sets of octahedral sites located at body-center locations (f ig. 1). It 
is related to the B2-NiAl structure in that every other Al site in 1 he N i A 1 lattice becomes 
occupied by a Ti atom in an ordered fashion. Because of the similarity in lattice structure 
between B 2 and L2 j . it is possible to develoj) ])recipitation hardened alloys sitnliar to /y ' 
nickel base superalloys. However, in order to design within this new family of B2/F2| 
alloys (also known as A/.T superalloys), microst ructural factors such as the solubility of 
Ti in NiAL the defect structure within the solid solution NiAl-Ti phase, the lattice misfit 
bdween the two phases, and eventually the effect of quaternary and higher order additions 
on all of these factors need to be determined in detail. However, even with the amount of 
experimental work performed on the ternary NiAl-Ti system over the last decade [2-4.25-2N], 
details of these microst ructural features are still vague while those for quaternary systems 
are completely lacking. 

Needless to say, theoretical modelling of these microst ructural features would save years 
of alloy development time and result in significant cost savings. In fact, in today’s climate 
of dwindling resources, the only feasible way to develop new alloy systems may be through 
advanced screening by theoretical means. However, almost no computational effort has 
been spent analyzing this issue. Fu et al. [29] have recently presented a comprehensive 
first principles study of site substitution on both the defect structure of FeAl and N i Al 
compounds as well as the behavior of ternary (Ti or Cr in FeAl and Fe in NiAl) additions 
to t hese systems. Unfortunately, no NiAl+Ti results were included. In a much earlier study. 
Tso and Sanchez [30] performed some thermodynamic modelling of the Ni-Al-Ti system but 
the results focused on the modelling of the various phase diagrams with limited attention 
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and no detail presented concerning the case of ternary alloys. 

In this paper, we concent rate on Ni-rich NiAl alloys, prevalent in industrial applications 
[3.4]. tor binary alloys, the defect structure consists primarily ol substitutional antisite 
defects on the A1 sublattice and as a result the lattice parameter monotonically decreases 
with increasing Xi concentration [1]. In a previous application of the HI S method to 
1 he study of the zero temperature defect structure of NiAl alloys [31]. 1 lie computational 
results correctly identified the energetically favored defect structures, including the correct 
substitutional defect scheme in Ni-rich alloys. It was also shown that while Ni vacancies are 
mostly responsible for the observed behavior in Al-rich alloys, evidence was found for the 
possibility of vacant sites in both sublattices as well as a Irend for clustering of vacancies. 
Moth features have been observed in recent experiments [32]. 

In particular, the quantitative accuracy of the BFS results for binary NiAl is highlighted 
by the predicted values of the lattice parameter. A survey [1] of available experimental data 
indicates that the lattice parameter in the Ni-rich region varies linearly with concentration. 
Correspondingly, the BFS results can be adjusted to a similar linear regime with a very 
small departure from the experimental expression: 





1.03X4 - 0.000791 4. rv, 
1.0390- 0.0007()X9.r\ ; 


(M) 


Vet) / b hs 

Fig. 2 displays both the experimental values as well as the two linear expressions pa- 
rameterizing both experimental and theoretical results. The qualitative and quantitative 
agreement with experiment certify the validity of the NFNi. A I - A I and Ni-AI BFS pa.ra.in- 
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eters used in this work. From a mot hodological point of view, the results that follow can 


bo taken as evidence for the validity of the Ti-Ti. Ti-AI and Ti-Ni BFS parameters used in 
this paper. 

A. Analytical Structures 

Fxploiting the computational simplicity of semiempirical methods and in particular 
the BFS method, this section is devoted to a discussion of the energetics of a large num- 
ber of NFAl-Ti alloys. These alloys represent a wide range of concentrations and different 
atomic distributions with different degrees and types of ordering. For reasons of simplicity 
and to enable comparison with available experimental data on industrially relevant com- 
positions. we restrict our calculations to Ni-rich Ni.^jAFo-./ri.r alloys. The large size of 
the computational cells combined with the number of different elements considered results 
in a large number of possible distributions of those atoms within the computational cell. 
Consequently, an efficient and economical - but physically sound - technique is required to 
examine all pertinent configurations of atoms. 

In this section, the approach used is half way between a detailed analytical study 
and an exhaustive statistical survey for each concentration and atomic distribution. The 
procedure followed was to define a particular ‘configuration* of atoms, then use the BFS 
mot hotl to compute the energy of formation of each configuration, as well as its equilibrium 
lattice parameter and bulk modulus. A large set of high symmet ry configurations covering 
a wide composition range were defined in this manner. The purpose of this approach is 
to l ) relate general trends of the lattice parameter and the energy of formation to changes 
in concentration and atomic distribution and 2) identify melastable structures, that is. 
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configurations with energy close to that of the ground stale, which might have a high 
probability of appearing in the actual alloy depending on tl ie conditions prevalent during 
its processing. Therefore, this ‘catalog' of alternative configurations serves the purpose 
of identifying ordering trends and provides physical insight to the microst ruct ure of the 
actual alloy. Since the number of possible configurations is prohibitively large, we believe 
that this selection procedure is sufficient to gain an adequate understanding of the alloy 
studied. Some selected configurations are described in Fig. T 1 he 72-atom computational 
cell is properly labeled to facilitate the understanding of the atomic exchanges that define 
the different configurations in the set (Fig. 5. a). Two examples of such configurations are 
shown in Fig. Tb. and a complete list, of the configurat ions studied in this work is included 
in Appendix B. 

The obvious difficulty in defining an appropriate set of configurations for the study at 
hand is the lack of a priori knowledge of the ground state structure for a given composition 
or temperature. In some cases (i.e. low concentration solutes) the ground state might 
have easily predictable features but. in general, no guarantee exists that the ground state 
configurations for each concentration are included. However, a close examination of the 
results usually give a clear indication whether a particular configuration could be 'missing 
from the set. 

In this work, the results for a set of over 150 different configurations of a 72-atom 
computational cell are shown and discussed. Fig. 5 shows some of these configurations 
and introduces the notation used in Appendix B to identify the position of Ni. A1 and Ti 
atoms in the lattice. Obviously, if this set. was complete and if it included the ground state 
for a certain composition, the energy of formation for such a configuration would be lower 
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than any other with the same concentration of elements. Fig. 4 summarizes the results 
of this type of analysis for the 150+ set of configurations: the energy of formation at zero 
temperal lire for NUoAUu-./I fr alloys, for .r between 0 ( B2 NiAl) and 50 ( B2 XiTi) at. 
*/<' I i. is shown in Fig. 4. a. A close examination of the atomic distributions considered 
for each concentration - detailed in Appendix B- reveals that most of the likely ordering 
patterns for a given composition are included in this survey. Moreover, Fig. 4 shows, for 
each concentration, ‘clusters' of states very close in energy that correspond, in terms of the 
atomic distributions, to related ordering patterns. In particular. Fig. l.b focuses on those 
configurations with 0(.r7 ,{*25. highlighting the fact that a specific group of configurations, 
corresponding to the lowest energy states for each concentration above xj; ~ 5 follows a 
trend quite different from the one that characterizes the rest of the configurations in 1 lie 
sot. In what follows, we first discuss the BFS predictions for preferential site occupancy of 
l i in NiAl. followed by an analysis of phase stability, based on the results summarized in 
Tig. 4. 

B. Site preference of Ti in NiAl alloys 

In general, those configurations where Ti atoms occupy sites on the A1 sublattice are 
consistently lower in energy. The site preference of Ti atoms can be analytically determined 
with the BFS method by considering just a few configurations where, for a fixed concen- 
tration, Ti atoms are located in specific sites. Figure 5 shows the configurations considered 
and Table 4 lists the corresponding values for the energy of formation and lattice parameter. 

The notation used to indicate site preference is self-explanatory and is convenient for 
representing more complicated situations as will be seen in an upcoming study of multicom- 
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ponent systems. The notation A(B) indicates a substitutional defect where atom A occupies 
a B site. The symbol A(H)(- indicates an A atom occupying a If site with the displaced B 
atom occupying a (' site. When necessary, the separation between 1 wo defects is noted : lor 
example. AtBJC.wy indicates that the A atom in the If site 1 and the B atom in the < 1 site are 
nearest neighbors. In Table •$. the subindex far (d) denotes that the corresponding atoms 
are at a distance d greater t han second neighbor distances, while t he subindex fl( usUr in- 
dicates that the atoms locate themselves in A1 sites following the Ileusler pattern (atoms 
in any pair of Heusler sites are located in opposite corners of a face of the elementary cube 
in the A1 sublattice). 

The cases considered for the study of sile substitution include two different Ti con- 
centrations: i.:f9 and 2.7K at. %. For the first case, two basic configurations are possible: 
the Ti atom occupying an Al site ( T i ( A 1 ) ) (fig. 5. a), and the Ti atom occupying a Ni 
site with the displaced Ni atom located in an Al site (Ti(Ni)AI). In this second case, the 
antistructure Ni atom can be a nearest neighbor of the Ti atom ( 1 i( Ni )Aly,v* fig- b.b ) 
or not ( Ti(Ni)Al/ fU . fig. 5.c). The corresponding values for the energy of formation are 
listed in Table T clearly showing the preference of a Ti atom for an Al site. A similar 
situation is observed at a higher Ti concentration (figs. 5.d-5.h). Several possibilities now 
exist concerning the relative location of the additional Ti atoms: they could both be lo- 
cated in Al sites, or create an antistructure Ni pair by occupying Ni sites with the Ni atoms 
occupying Al sites, or a combination of these defects with several options for the relative 
location of the antistructure and substitutional atoms with regard to eacli other. All these 
possibilities are included in Fig. 5 and the corresponding formation energies are listed in 
Table d. Once again there is a clear preference of Ti for Al sites compared to any other 
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substitution scheme. Moreover, the small difference in energy of formation between the case 


when the two Ti atoms are nearest neighbors and the case whore they are truly separated 
can be used to explain why at low Ti concentrations there is a close competition between 
a disordered and ordered arrangement of Ti atoms on the Al-sublattice. as seen in Fig. 
l.b (ordered configurations are indicated with circles while disordered configurations are 
indicated with solid squares). The preference of Ti for Al sites is observed for the whole 
range of Ti concentrations considered in this work, leading eventually to the formation of 
the Heusler ( N i > AlTi ) alloy. 

Figure (i is an alternative way of displaying the results of the site preference calculations, 
in the form of an energy spectrum. The first column shows the energy ‘levels' corresponding 
to the different substitution schemes for just one Ti atom in a. 72 atom cell. The second and 
third columns show results for two (,rj, — 2.78) and five ( xj t — 0.94) atoms respectively. 

The splitting of the Ti( AI) level for 7 = 2.78 is due to the different relative locations of 

the two Ti atoms: they could be 11011-interacting (‘far'), situated in a Heusler-like pattern 
(opposite corners of a cube in the Al sublattice, indicated wit h circles in Fig. 4.b) or in such 
a way that the two Ti atoms in the Al sublattice are at next -nearest -neighbor distance. For 
x/v =(>.94 the splitting in energy levels is much more since many more possibilities exist 
for the placement of the Ti atoms. We only show two states in the ground state region: 
the one corresponding to the Ti atoms in solid solution ('far' from each other in the Al 

sublattice) and the Ti atoms following a strict Heusler ordering pattern. The reversal 

in energy levels between X77 =2.78 and X77 =0.94 for the solid solution and the Heusler 
ordering options clearly indicates that the formation of Heusler precipitates is favored at 
t he higher Ti concentrations. 
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X I i 

Description of 

Energy of 

Lattice 


the configuration 

Formation (eV) 

Parameter ( A ) 


Ti( Ni )AI\\ 

-0.4 1420 

2.X 70 

1 .49 

Ti( Ni )Al/„, 

-(U,sor>2 

2.S70 


Ti(AI) 

-0.()0X.49 

2.X 54 


2Ti(.\i),v,v + 2Ni( Al) v \ 

-0.:{20()7 

2. XX 1 


2[77( .V/}.l/].v v 

-0. 42X74 

2.XX2 


Ti(Ni)Al v v + Ti(Al) 

-0.4()X94 

2.X 79 

2.7X 

[ / i( A i 1 + A /(.!/) + / i{ . 1 / )]y 

-0.47774 

2.X7X 


2Ti(Al)v.v 

-0.599 IX 

2.X59 


2Ii( A 1 )j ar( y y2 ll j 

-0. 00450 

2.X5X 


2Ti(Al)/y f 

-0.00450 

2.X5X 


2J i(Al) /(ir(v yr{„) 

-0.00404 

2.X5X 


Table 4: Energy of formation and lattice parameter for specific atomic configurations of Ti 
atoms within a NiAl lattice at two different concentration levels. 

Finally, the stability of the Heusler phase (Fig. 4,b). where Ti atoms are located 
exclusively in A1 sites, indicates that the site preference observed at low concentrations is 
the same for t he whole range of concentrations for which t his phase exists. 

C. Evolution of second-phase structures within NiAl 

Figure 4 introduces some selected configurations from the complete set used in this 
work, as well as the notation used in labeling them. A complete list of configurations 
is included in Appendix R and the energy of formation computed via BI S for all the 
configurations is shown in Fig. 1. One feature in Fig. 1 becomes immediately apparent: 
beyond 5 at,. ( /l Ti. there is a clear separation in formation energy between a selected type of 
configurations (denoted with circles in Fig. 4.1)) and the rest (denoted with solid squares). 
The selected group of configurations correspond to a particular type of ordering, where Ti 
atoms locate themselves exclusively in A I sites in such a way that they always have Al 
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atoms as next nearest neighbors. This ordering pattern, when extended to the case when 
the concentration of Ti is *25 %, corresponds to the L2i ordered structure or Heusler phase. 
At that particular concent rat ion. the energy gap between the Heusler structure ( solid circle 
in Fig. Lb) and any other configuration is the largest, clearly indicating the stability of 
this phase at a stoichiometry equivalent to Ni-^AITi. 

Melow 5 at. ( Z Ti. those configurations where Ti is in solid solution with the matrix 
are energetically favored, however so slightly, over those where short range order dominates. 
This is more clearly demonstrated in Figure (j. This situation is reversed at a composition 
near 5 at. l Z Ti. with an ever more distinguishable preference for Heusler ordering against 
any other option beyond 10 at. { Z Ti. In other words. Heusler-like ordering becomes clearly 
preferred beyond a certain critical value somewhere near 5 at. ( /f Ti. in spite of the fact that 
even at lower concentrations those configurations with short-range Heusler-type ordering 
an' also very low in energy. This crude way of determining the solubility limit of Ti in NiAL 
which based on these results could be set somewhere near 5 at. % Ti. not only establishes 
such a critical value but it also provides some insight on the behavior of the system for a 
wide range of concentrations surrounding the solubility limit. 

One clear reason for the stability of the Heusler phase is the typo of nearest neighbor 
bonds present. The Heusler ordering maximizes the number of energetically favorable Ni-Al 
and Ni-Ti bonds - both B2 compounds - as seen in Table 4. In addition, the second neighbor 
bonds are either Ti- Al. also energetically favorable, or Ni-Ni. The \i-Ni bonds, due to the 
closeness between the lattice parameter of the Heusler phase to the equilibrium value of the 
lattice parameter of Ni in its bcc phase, introduce very little strain in the lattice. 

A closer examination of the configurations in Appendix B shows that the preference 
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for 1/2 1 ordering is also apparent in the high energy states in the nonseleeted configurations 
in Fig. |. This is because any configuration that contains Ti atoms in Al sites wit h only Al 
next nearest neighbors introduces a significant gain in energy. An additional advantage of 
this type of analysis is that a good deal of information on alternative st ruct ures can also be 
obtained. This is an important issue when dealing with metastable and other alternative 
structures which may result from a particular processing scheme. For instance, we could 
examine the energetics of an alternative type of ordering, which for example could involve 
a slight change in site occupancy in the Al-sublat tice. Fig. 7 shows the Heusler unit cell as 
well as one corresponding to this alternative ordering scheme (to be called Kneen phase), 
which shares the first neighbor coordination wit h the Heusler phase. The only difference in 
these structures resides in the second neighbor coordination, with two I i-Al bonds being 
replaced by pure Ti-Ti and AFA1 bonds, due to the different distribution of Ti and Al 
atoms in the Al-sublat tice. Due to the size of Ti and A! atoms relative to the lattice 
parameter of the Heusler phase, this change introduces enough strain in the lattice to make 
this configuration energetically less favorable and therefore less likely to be found. However, 
in a realistic situation and depending on the processing conditions of the alloy, it would not 
be improbable to find short range order similar to that of the Kneen phase (Fig. T.b) in 
addition to the expected Heusler ordering. Moreover, in regions with low Ti concentrations, 
the Kneen phase competes with the Heusler phase in that they both share a great deal of 
common short-range order. 

The somewhat large difference in energy between configurations with the same con- 
centration is not reflected in the corresponding values of the lattice parameters (Appendix 
B). The lattice parameters of the different ternary alloys considered follow a nearly linear 
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Atom 

Nea 

rest 

neighbors (a), (b) 

Next - nearest -neighbors ( a) 

Next -nearest - neighbors ( b ) 




li 

Ni 

Al 

Ti 

Ni 

Al 

Ti 

Ni 

0 

4 

4 

G 

0 

0 

G 

0 

0 

Al 

B3 

0 

0 

0 

0 

G 

0 

2 

4 

Ti 

IB 

0 

0 

0 

() 

0 

0 

4 

2 


Table I: Coordination of nearest -neighbor and next-nearest -neighbor bonds for Ni. Al and 
Ti atoms (lirsl column) in the (a) Heusler phase and in (I)) Kneen phase. as shown in Fig. 

i . 

relationship between the binary Ni Al and Ni li B2 values as a function of Ti concentration 
(Fig. N.a). There is. however, a small deviation below the ideal linear relationship, common 
to most of the configurations that display short- and long-range order. The overall trend, 
however, is clear. The lattice parameter of ternary NisoAl^o-.rTi.,. alloys increases with 
increasing Ti content, as confirmed experimentally in Ref. 4. A more detailed display of 
t hese results is seen in Fig. N.b, where the lattice parameter for the configurations in Fig. 
4.1) are shown. 

Table 5 lists the BI S predictions for the lattice parameter of the lowest energy config- 
urations for each concentration. Two values are shown: 'Ordered', indicating the configu- 
ration where Ti atoms locate themselves exclusively in Al sites following Heusler ordering, 
and 'Disordered', where the Ti atoms are randomly scattered in the Al sublattice. 

The re suits of the first column ('Ordered') can be easily adjusted to an expression of 
the form 

— — = 1.00000 + 0.09:M7.r r , 0 < x Ti < 25 ( F r >) 

<*X,Al 

and is plotted in Fig. N.b, indicating the linear dependence of the lattice parameter of 
NiAl-Ti alloys with increasing Ti concentration. 


NASA TM-1 13121 


22 








.r/-, 

Ordered 

Disordered 

2.7X 

2.X50 

2.X50 

4.17 

2.X59 

2.X59 

5.55 

2.X03 

2.<s(j:$ 

(j.94 

2.X07 

2.XG7 

X.33 

2.X71 

2.X 70 

9.72 

2.X 75 

2.X74 

11.11 

2.X 79 

2.X 7 7 

12.50 

2.X S3 

2.XX1 


Table 5: Lattice parameter (in A) of the lowest energy NiAl-Ti configuration for several 
li concentrations (see Appendix B (or lattice parameter values for the complete set of 
configurations used in this work. 


Alternatively, the BF rule [19] - a simple rule derived from the BFS equations - can pro- 
vide an expression for predicting the lattice parameter for NiAl-Ti alloys ol any composition. 
Using Eq. (39) of Ref. 19. and the parameters listed in Table 1. we obtain 




lXX7.32.r iV ,+ 79U3.r, w T 1249.44.r T ; 


(lb) 


(jX5.X0.rv; T 24X.XX.r 4 / T 3 XX.s 7 .r 7 1 , 

Similar expressions (Eqs. (40)-(47) in Ref. 19) can be obtained for the concentration 
dependence of the bulk modulus and cohesive energy per atom. Also, simple rules can be 
obtained (Eqs. ( 45 )-( 47 ) in Ref. 19) regarding the deviation of such values with respect to 
the usually assumed average values. 


D. Monte Carlo Simulations 

While the information provided by st udying a large number of candidate configurations 
at absolute zero temperature provides valuable information on the energetics of the system 
at hand, it says very little in terms of the effect of temperature and processing on the 
microstructure of the alloy. In this sense, several numerical tools exist to investigate these 
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issues. In t his work, wo concentrate on t he application of Monte ( 'arlo procedures to perform 
simulations of Ni-Al-Ti alloys. When added to the results of the previous section, the 
simulation results should provide a more complete theoretical pict ure of the microst met ural 
behavior of these alloys. 

The Monte (’arlo procedure employed in this work is a variant of that used to simulate 
the Ising system. The computational cell used in all simulations consists of a large number 
of atoms (1021 in most cases shown in this paper) arranged on a body-centcred-ciibic lat- 
tice. Boundary effects are minimized through the use of periodic boundary conditions in all 
directions. Although the simulations shown correspond to a wide range of temperatures, 
for simplicity we ignore the lattice parameter dependence on temperature and use the T= 
0 predictions discussed in the previous section, as we are mostly interested in ground state 
the computational cell is a random solution, i.e. the Ni. A1 and Ti atoms are randomly 
assigned to each site, in direct proportion to their assigned composition (unless otherwise 
indicated). No vacancies were allowed in the majority of simulations reported in this work, 
though their incorporation does not present any added degree of difficulty, as will be demon- 
strated shortly. Most of the simulations shown correspond to a process where a sequence 
of decreasing temperatures (the ‘cascade*) is chosen, where the system is allowed to equi- 
librate sequentially at each temperature. This simulates the ‘slow cooling* of the actual 
alloy, starting from a high temperature disordered solid solution. The equilibration pro- 
cedure involves the random selection of a pair of atoms and their subsequent, reversal in 
position. The reversal in chemical identity is accepted or rejected using the Metropolis crite- 
rion where, for a given temperature, the probability assigned to a particular exchange A — 
B is ( Ah i kl , where k is Boltzmann constant and A A is the change in BPS energy between 


NASA TM-1 13121 


24 



1 1 io configurations after and before t lie switch. I he number ol switches allowed is set as an 
input parameter and it determines, toget her wit h the difference in temperature between two 
successive steps in the simulation, the simulated cooling rate of the sample. A large number 
of switches is equivalent, to a longer stabilization time, which together with small changes 
in temperatures simulate a slow cooling rate. In this calculation wo ignored the dependence 
of the lattice constant with temperature, therefore temperature enters into the calculation 
only through the Metropolis criterion. After the system has achieved equilibrium (based on 
the total energy of the computational cell), various properties of the system are computed 
and averaged and then the simulation proceeds to the next temperature decrement. The 
properties calculated include the average energy of the cell, the specific heal, and bond 
correlations. 

While these simulations do not attempt to mimic the detailed dynamics of the equili- 
bration process, they do offer a qualitative view of the effects of rapid versus slow cooling 
of the system. The cooling rate ( that is. the size of the steps between the various tem- 
peratures considered in the cascade) is of critical importance in determining the final state 
of the system. Slow cooling results in a highly-ordered low-temperature state, while rapid 
cooling results in a more disordered material often containing antiphase boundaries. As 
with actual processing, the temperature treatment of the sample is essential in determining 
the final state. It is to be expected that the slow cascade processes, used in most of our 
calculations, will result in highly ordered compounds. In contrast, sudden cooling of the 
sample will result in regions of disorder often in the form of antiphase boundaries . It is also 
possible that additional phases - like the ones described in t he previous section - might also 
appear with a frequency proportional to the difference in energy of formation with respect 
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toother structures. The presence of these higher energy structures in the fast cooled sample 
can arise if quenching of the sample "freezes' specific domains within the cell, whoso seed is 
already present in the initial disordered state. 

An example of the effect of cooling rate on structure is shown in Fig. 9 which includes 
the results of two separate simulations for an alloy with a bulk composition equivalent to 
Nij Al i i, both starting from the same random state at high temperature. To facilitate the 
visual interpretation of the results, the 1024 cubic cell is stretched along the 001 direction. 
The rapidly cooled cell, shown in Fig. 9. a. is characterized bv order and disorder features: 
there is a clear trend towards separation of Ni-rich and Ti-Al alternating planes, as well 
as some indication of Heusler ordering in about half the sample. The lower half of the 
cell shows a noticeable degree of disorder and even the existence of some energetically 
unfavorable situations, mostly in regions of large Ti or A1 concentrations. Figure 9.1). in 
contrast, shows the results of a cascade process that had finer steps between temperatures 
(slow cooled). Kven though the final temperature was the same as in the previous case, this 
cell displays almost perfect Heusler ordering. Most of the antisite defects present in this 
sample are eliminated after subsequent re-heating and slow re-cooling - a second cascade - 
which results in an almost ideal Heusler phase with a few antisite atoms. 

Figure 10 displays similar results for the slow cooling of a Nir l0 ( ALTi)no alloy with 1 at. 
{ A Ti. Starting with a random alloy at high temperature, the slowly cooled sample settles 
into a perfect H2 NiAl ordering with Ti atoms in Al sites. However, no trend toward pre- 
cipitate formation is seen at this concentration, consistent with the configurational analysis 
described in the previous section. 

Figure 11 shows results for a Nir l0 ( Al.Ti )r>o alloy with 5 at . V( Ti for two different 
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computational cells. The first three cells in Fig. 1 1 show the final results of three consecutive 
temperature cascades, each showing a higher degree of short range order. 1 he last cell 
displays 1 he results of a single temperature cascade for a large ( lf>()N atoms) cell, showing 
essentially the same features observed in the smaller cell. Since this concentration is near 
the accepted solubility limit for Ti in NiAl, a larger cell allows for a closer examination 
of the final state and the possibility of precipitate formation. In this case, t he Ti atoms 
actually display three different types of behavior, as seen in the last column in Fig. 1 1. which 
corresponds to the final temperature (T = 100 K): 1 ) the format ion of Heusler precipitates, 
more clearly seen if the periodicity of the cell is taken into account . 2) Ti atoms in solid 
solution in the NiAl matrix with preference for the Al-sublal t ice. and 4) the distribution of 
Ti atoms following the Kneen ordering (atoms located in alternate sites along rows parallel 
to the edges of the sample ) in regions often adjacent to Heusler precipitates. Also, a few 
Ni antisite atoms are seen (red atoms in the blue planes). While it might prove to be a 
premature conclusion, it is noticeable that, in the final state the antistructure Ni atoms seem 
to 'attract" Ti atoms to their vicinity in the A1 plane, creating a Ti- rich interphase with the 
NiAl matrix. 

Similar struct ures are observed in a. Ni-40Al-10Ti alloy which has undergone a double- 
cascade process (Fig. 12). At the end of each cascade, the formation of Heusler precipitates 
(characterized by the chains of Ti atoms along diagonals in the A1 planes) is clearly observed. 
In spite of the re heating and re-cooling of the sample, the Heusler precipitates reappear 
proving that the dominant effect of sufficient additions of Ti to Ni-rich NiAl is the formation 
of such precipitates. Some of the other features pointed out in Fig. 11 are also featured 
in the 10 at. % Ti case: the solid solution of Ti atoms in the matrix, and the presence of 
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anl istructure Ni atoms in the A1 sublaltice, with an apparent clustering of Ti atoms in the 


vicinity. Also, some Ti atoms following the Kneen ordering are clearly seen. 

f inally, we would like to comment on the role of vacancies on the microst met lire of the 
alloys studied. So far. all the examples shown ignore the presence of vacancies. If vacant 
sites are allowed in the calculation, it is found that none of the essential microstructural 
lea til res discussed earlier are affected to any significant degree. However, the results show 
clear indication of vacancy clustering in such a way t hat inner Al surfaces are created. This is 
not surprising, given the large size of Al atoms and the low surface energy of Al, both features 
conducive to the formation of Al surfaces. Moreover, due to the tendency of vacancies to 
coalesce. Ni ant istruct ure atoms are t hen found in the vicinity of vacancy clusters. Some 
of the unfavorable Ni-Ni bonds thus created are compensated by the migration of some Ti 
atoms to Ni sites creating favorable Ti- Al bonds. 

IV. Experimental Analysis of NiAl-Ti Alloys 

Three Ni Al single crystal alloys (Ni-47Al-3Ti. Ni-45Al-5Ti and Ni-43Al-7Ti) were 
grown by a Bridgman technique at the University of Florida. The ingots were homoge- 
nized for 32 Ill’s, at 1(>44 K, aged for (> hrs. at 1255 K. and slowly furnace cooled from the 
aging temperature. The purpose of this heat treatment was to produce a low temperature 
'equilibrium' microstructure that would best correspond to the ground state condit ions mod- 
elled under the BPS technique. Samples for transmission electron microscopy (TEM) were 
prepared from 3 mm diameter cylinders electro-discharge machined from the heat treated 
ingots. Slices sectioned from the cylinders were mechanically ground and elect rochemi- 
cally t hinned in a twin-jet Tenupol 3 polisher. Microstructural ( bright -field /dark-field ) and 
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diffraction analysis wore conducted in a Phillips 1001 I KM equipped with a double till 
goniometer. 

Fig. Pi. a shows a hrighl field image of t he X i- 1 7 A 1-3 l i alloy. Except for an occasional 
dislocation, t he microstructure is very clean and featureless and shows no sign of any second 
phase precipitation. 'Phis is confirmed in the corresponding { 1 JO) zone-axis selected area 
diffraction pattern (SAI)P), presented in 1 ig. Pi. I k which shows only the NiAl matrix spots 
and no extra diffraction features (spots or streaking) due to precipitation. 

In contrast. Fig. 14. a shows a bright-field image of the Ni-45AI-5Ti alloy after the 
same thermal treatment. Precipitation of a high density of extremely fine second-phase 
particles can be clearly seen, especially in the dark-field image shown in Fig. I F!), where 
the precipitates appear bright on a black background. The corresponding (110) SADP in 
Fig. l l.c shows distinct extra spots which were indexed to a fee crystal structure with 
lattice parameter a 0 = 5.<S(j A. corresponding to the Ni 2 \ITi (Heiisler) phase. From the 
crystallographic information revealed from the diffraction pattern (schematically illustrated 
and labeled in Fig. 14.d), it can be seen that the Ni 2 AlTi phase nucleates with a cube- 
on-cube orienlalion relationship with the NiAl matrix, i.e. [1 10]/y,4///[1 1()]ak 2 .i/t< and 
(001),v;.4///(001)a’/.>.4/T4- Due to the fine size of the precipitates, 1-5 mil, and the small 
lattice misfit between these two phases, on the order of 1.5 ( X , the precipitates are coherent 
with the matrix resulting in significant coherency strains around the particles, as seen in 
t he bright-field image in Fig. 1 Fa. 

Figure 15 shows a dark-field TEM image of the Ni-43A1-/Ii alloy after the same heat 
treatment. Precipitation of a high density of well defined and coherent Ileusler precipitate 
plates ranging in size between 10-50 mil can be clearly seen. While tin 1 precipitates are on 
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average 10 times larger than those in the \ T i-45AI-- r )Ti alloy, due to the small lattice misfit 
between the j)reci[)itate ph asi' and the NiAI. most of t lie plates are still coherent with the 
matrix. 

V. Discussion 

In exploiting the computational simplicity of the MI'S method in calculating the en- 
ergetics of various alloy structures, it must be understood that the technique does not 
automatically provide the ground state configuration of a particular alloy but depends on 
the inclusion of that configuration in the catalogue of configurations selected for study. 
For reference in this and future studies, the present catalogue of configurations has been 
defined in Appendix IF Returning one last, time to the issue of completeness of this set of 
configurations. Fig. 1 can be used as a clear example of how an occasional omission in the 
predetermined catalogue can be easily detected and corrected as the large size of the set 
provides enough statistical information to detect trends and patterns and therefore, infer 
what configurations might be missing. Had one of the energetically favorable configurat ions 
been left out of the original set, the omission would have been noticed as a "discontinuity ' 
in the plot shown in Fig. 4.b. For example, had the Ni-jAlTi Heusler phase (denoted with 
a circle in Fig. 44>) not been included in the set. it would have been easily inferred by 
observing the series of states indicated with open circles, leading to the prediction of such 
a phase as the ground state for that concentration. In fact, the analytical MFS results for 
the 25 at. % Ti alloy not only suggest but confirm the I/2i structure as the ground state 
for Ni.AlTi. 

The use of this survey met hod for studying t he energetics of various alloy systems pro- 
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vidos one wit-li significant information on not only the ground state structure oi a particular 
alloy but relates general trends in the energy of formation and lattice parameter to changes 
in concentration and atomic distribution. It also has t he important additional advantage of 
being able to identify metastable structures or configurations with energy close to that of 
the ground state, which may have a high probability of appearing in the alloy depending on 
the actual processing conditions. We note also that Monte Carlo simulations complement 
the use of our survey method to obtain ground state structures. 

One of the more significant results of this investigation is the successful application of 
the BFS method to the determination of the solubility level of a ternary addition to an 
ordered intermetallic compound, in this case Ti in NiAl. The results of static calculations, 
shown in Fig. 4. suggested that the solubility of Ti in NiAl is ~ 5 at. %. Ibis was 
experimentally verified by examining the microstructure of three NiAl-Ti alloys, where 
it was found that Ni-47AI-3Ti was a complete solid-solution alloy with no second-phase 
precipitation, and the nucleation of a high density of N^AH i precipitates appeared only in 
NiAl alloys containing 5 or more at. % Ti. Although the calculation of minimum energy 
ground state configurations via the BFS method was performed at 0 K and t he experimental 
results were obtained at room temperature, the microstructure of these high-melting point 
ordered alloys is not expected to be too different at these low temperatures. 

Another useful result of the BFS calculations is the ability of the method to predict 
the atomic structure of the second phase particles and the lattice parameters of all the 
constituent phases. Both the static and Monte Carlo simulations were able to predict the 
correct atomic configuration of the resulting second phase particles, i.e. N^AlTi (Heusler 
phase), in NiAl-Ti alloys. Based on these calculations, the lattice parameter of the Heusler 
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phase was found 1o be 0.5828 did, which is in close agreement with the average value of 
0.5X70 mu reported for the NFAlTi st ructure [55]. The lattice parameter of the correspond- 
ing solid solution NiAl-Ti alloy was calculated to be 0.2X05 mil. Tims, based solely on the 
BIS analytical approach, the lattice mismatch between N i A I and N i j AIT i is found to be 1.7 
( /<A which is very close to the approximately 1.5 % misfit obtained experimentally via TEM. 
A controlled mismatch between phases is one of the major design criteria in almost all high 
temperature alloys. Therefore, the ability to model not only the correct second phases but 
also the resulting lattice mismatch is a significant breakthrough in the computational design 
of high temperature alloys. 

An additional significant contribution of the BFS method is the ease with which the 
site occupancy of a third element can be determined in a structure. Knowledge of site 
occupancy of an alloying element is an important and necessary piece of information in 
understanding the defect structure and its impact on mechanical properties [54]. In the 
case of NiAl-Ti alloys studied here, the BFS method was able to correctly predict the site 
occupancy of Ti in the NiAl lattice. The preference of a Ti atom to occupy an A1 site in 
Ni rich NiAl is well established [4] and. in fact, all the empirical alloy development programs 
start with this basic assumption. 

The most common experimental techniques for determination of site occupancy are 
A hr HEM I (atom location by channeling enhanced microanalysis) [55] and AFFIM (atom- 
probe field ion microscopy) [5(>]. However, these are often tedious and very involved pro- 
cedures and in the case of the former there are a number of complications which make the 
technique relatively inaccurate [57]. The ability and ease to analytically determine the site 
occupancy of alloying additions, as proven by the BP'S simulations shown in this work, is 
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not only beneficial to an alloy developer, but is a boon to those trying to determine the 
mechanistic behavior of ordered alloys. Moreover, the methodology used in this work is 
easily applicable to more than one alloying element , which allows for a clear understanding 
of the interaction between the different alloying additions. 

While a great deal of structural information has been derived lrom Ihe static calcu- 
lations. they do not provide much information in regards to the effect of temperature and 
processing conditions on microstructure, therefore, to complement the static calculations 
and furt her our understanding of these alloys. Monte Carlo procedures were used to perform 
simulations on the Ni-Al-Ti alloys. These results are summarized in Fig. 10. Figure 10 
shows the 1024-atom computational cell for the final states of temperature cascades for a 
number of different concentrations, some of these have been discussed in detail previously. 
The results show that Heusler formation is apparent beyond the solubility limit for Ti. But 
as in anv real system, there is also a statistical chance for t he development of other struc- 
tures at a rate based on differences in energy between them and the ground state. In this 
case. Kneen ordering and an occasional Ni-antistructure atom are evident . The abundance 
of Kneen ordering can be easily explained in terms of our earlier discussion ol nearest and 
next-nearest -neighbor coordination. With such small energy differences involved, when a 
Ti atom has the choice of occupying an A1 site in a Monte Carlo calculation, the difference 
between a Heusler site and a Kneen site makes the probability at high temperatures almost 
identical between the two sites. Therefore, once a Ti atom is "trapped' in a. Kneen site, its 
likelihood to migrate further to form a Heusler arrangement with decreasing temperature 
diminishes. With increasing Ti concentration, the energetically favorable Heusler site be- 
comes more prevalent and the formation of Heusler precipitates is clearly favored due to 
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t he larger number of energetically favorable bonds t hus created. In addition, t he number of 
temperature decrements within the cascade - the cooling rate of the alloy - influences the 
density of these higher energy structures observed in the final microstructure. 

VI. Conclusions 

The BFS method has been successful in predicting the solubility limit, structure of 
the second phase particles, lattice mismatch between the alloy matrix and the precipitating 
phases, and preferred site occupancy of the alloying additions. These results have proven 
that almost all 1 he necessary parameters needed for a purely analytical alloy design approach 
are now within reach. The present results provide confidence in the BFS technique, the 
authenticity of the input parameters used (Tables 1 and 2) and the approach used for 
obtaining parameters by use of ab-initio methods, which removes the limitations imposed 
by the otherwise required experimental data base. The results also provide an energetic 
description of the detailed microstructure of the NiAl-Ti system. The real asset of the BFS 
method would be to model more complex systems containing two or more alloying additions, 
required for modelling alloys needed in practical applications, which is in progress. 
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Appendix A - BFS Parameterization 


An advantage of the BFS method is that, duo to its simple formulation, it allows for a 
straightforward (analytical) determination of the BFS S\h and parameters, therefore 

avoiding numerical uncertainties inherent to any numerical fitting procedure. Moreover, the 
input data used, whether it is obtained experimentally or from other theoretical calculations, 
'localizes* the accuracy of the ensuing BFS predictions for those alloys in the vicinity of the 
phase diagram of the ordered structure used as input. 

For simplicity, we reduce the following derivation to the case where the ordered struc- 
ture corresponds to a cubic lattice characterized by a single lattice parameter a,, (simple 
bcc or fee alloys with no tetragonal distortion). 

Consider an alloy A B. where due to the symmetry of the structure, there are 
non-equivalent atoms of species X (X = A. B) and n\ t denotes the multiplicity of the / th 
nonequivalent atom of species X. so that 


IE "v, = v ( i) 

A i = l 

where N ( is the total number of atoms in the cell. In this case, the t wo conditions used to 
determine the BFS parameters consist of exactly reproducing the heal of formation A// 0 
of the ordered structure and the corresponding lattice parameter, which are determined via 
LMTO calculations for a given structure. In this work, we used the B *2 N i A 1 base alloy as 
the basis for the LMTO calculation. 


V / = 1 ‘ 1 


( 2 ) 
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and 




EE 

A* /=! 


'ft X, <)< X, 
Oa 


= 0 


( :») 


where ( v ( represents the BFS cunt ribut -ion to t he energy of formal ion. given by 


r Y, - f \ t + a, ~ (( \ t ) 


( 1) 


Iii K<). (A 1). f y is the HFS strain energy and s y, — f( \ t ~ (< \ l ' s clieinical energy. 

< ' v u being the reference energy. 

For bulk ordered alloys, such as B2 NiAl. ihe HIS si rain energy is the same tor all 
atoms of the same species and it is uniquely determined by the input value of the lattice 
parameter «-o. 

f\, = I 'r {!-(!+ «a’ )<J \, } ( r >) 


where 

«A* = /t\'(«0 - «, A ) ( ( ») 

where and a} are the cohesive energy and equilibrium lattice parameter respectively 
for atoms of species X and — ({/l\- where q is a struct ure constant and l\ is a scaling 
length for species X [9.21]. The *glue* term included in Kq. (A 1) is given by 


tf.v, = ( 7 ) 

The HFS energy depends on the lattice parameter of the alloy structure only via the 
HFS strain energy and glue, therefore, Eqs. (A2)-(A3) can be written as 


E E ~rr ^ v ( ° u * + (JX ' ( H ° i A, } - ±n o ( x ) 

v ; 1 * c 
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"0 


v - ' y' ».v, i V, 

\ 2-* A',. 1 <)« 



<><lx, ( » ) 

<la 


} = 0 


Kt|. (AS) and (AO) can (lion lx* written as 


( 0 ) 


E E ^.v>o )^v. = a//„ - E E 


( 10 ) 


EE 0)c C Y, = EE A v « a* ( «o )U.x,( «o ) (11) 

V i x i 2 ' c 

If we concentrate now only on binary alloys that form fee or bee ordered structures 
characterized by a single lattice parameter ( L 1 2 - bio- B2, B32. etc.), then // 4 + //# — Ay. 
Simple expressions can then be obtained for the BFS chemical energies ? ( \ and 5^: 


and 


»Afl { A ] (‘ih ~ ;i A ) 

h ~ iUf >1 

- A. 4 ) 


wlioro 


= (" 4 + «h)-A// 0 - (W4< - 4(«o) + Wfl'fl(«-o)) 


( 12 ) 


( 13) 


( M) 


and 

h = E , 'V J A i 'rfl < v ) («a) ( 15) 

.V 

with try y )} — (}\ (do). The BFS chemical energies can then be determined with Eqs. (A 12) and 
(A 13). so that we can then search for the set of parameters ( A 4#, A# 4 ) that simultaneously 
satisfy these conditions. This is done by starting with the use of the following expression 
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for t ho BPS chemical energy in terms of the equivalent (chemical) lattice parameter 

S ( x = ' V I ‘-\ '! 1 - ( I + « v" I ( ( IB) 

where -) = 1 if <{ ( x ~ > 0 and “) = — 1 otherwise. The scaled lattice parameter for the chemical 
energy fly*, given l>v 

(i [x* = Av(«a- - ) ( 1 <") 

is rolatod 1« I ho HI’S parameters (_Yt h- ) •>>' moans of I ho HI’S equation for I ho chetniral 
energy 


+ Mrt'l x ( (ox + T 7 )/ ^ 


= E A ***3p“ (ox+A * x,rx ' + E '/vA'-y 


-(a v+^/>.V+ rrlcx') 


a x ,A j ( IS) 

A’ A- 

where N(M) is the number of nearest -neighbors (next-nearesl-neiglibors) in the equivalent 
crystal of species X. Y?i — c<7 y , i ?2 = fly , /'i — ca} and r\ 2 — (c = ^ for bcc). All 
three terms in the l.h.s. of Eq. (A IS) are known. 

Q v = V { o v i + + + A/R v + ^' X v+^)K A + ^-» ( 19) 

y\ o.v 


</i V = v,v ‘ 


20 ) 


and 


<12 = 


„ — ('■> A + XT“b 


( * 21 ) 


Because of the typical magnit ude of the ex])onent in Eq.( A21 ). it is reasonable to make the 
approximation ry 7 — r\ l only in that term, so that the l.h.s. of the BPS equation (A IS) 
reads 


Qx 


<f\ - f l 2 


= .\ XY r p x x f~ {t ' x+Ayx 


'■*1 


+ U.\ Yr 1 ?/ 


v x + v v h'-Vj 


( 2 * 2 ) 
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( **) 




This result is exact if second- neighbor contributions are not taken into account, or if for 
any particular structure M\y - 0 for all X and Y. as is the case in this work, where the 
B2 structure is such that any given atom has an atom of its own species as a next -nearest - 
neighbor, i.e. M y.v = 0 and My\ = 0. With the exception of the numerical solution of 
Eq. (Alb), the procedure for the determination of Ay-.v using Eq. ( A 2 b ) is straightforward 
and simple. Moreover, it can be easily shown that in most cases a quadratic approximation 
to the Rydberg function ( 1 + z)f~~ suffices to guarantee accuracy up tp 10 % of the exact 
results, with the added advantage of a completely analytical determination of the BFS 
parameters A and A^ 4 . The parameters used in this work were obtained by following 
the procedure described in this Appendix, including second neighbor interactions and a 
numerical solution of the trascen dental equations involved. 
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Appendix B - Atomic configurations 


The computational cell is defined in Fig. 2. a. It corresponds to a her lattice with 
7 2 sites. The B2 NiAl alloy corresponds to the atomic distribution shown in Tig. 2. a. 
where Ni atoms arc 1 denoted by black disks (labeled 1. 2. 2. ...) and the Al atoms are 
denoted by grey disks (labeled 12.14,12,...). A set of configurations is defined by changing 
the occupancy of these sites by 1) exchanging an A atom in site // with a B atom in site 
m ( A n — B m ) or by 2) substituting an atom B in site m with an atom A that originally 
was in site v (A n — li w ). Some of the configurations correspond to smaller versions of t he 

72 atom cell: those denoted with an * correspond to a cell where atoms In (n — 1 IS) 

have been eliminated and those configurations denoted with ** correspond to a fraction 
of the original cell where only atoms ( 1 .2,5.().9,1(M 7. IS. 25,2().29,20,27,2tt, 1 1. 12) are taken 
into account. Fig. 2.1) shows two examples of the configurations included in the set listed 
below, corresponding to .r/v = 11.11 and .rn — I C>.C>7. 
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1.17 

( . 1/ — / i ) 173). 71 

-0.6025X 

2.X59 


( .1/ — I i}\ 8.:t8.43 

0.00195 

2.X 59 


( 1/ - TiU 2-41 

*0.59592 

2.8(50 


( . 1 / — 1 / >42.48,47 

-0.59573 

2.XO0 


[Al — / / >38.42.48 

-0.59207 

2.XO0 


( A / — I / >5u.r>l.. r >8 4- ( *4/ — A / >82.88.70 

-0.41575 

2.X 73 

5.55 

(A/ / / ) 1 7, .47. 42 

-0.59X01 

2.X03 


( *4/ ~ ? ^37.42.47 

-0.59X44 

2.8(53 


( . 1/ — / / ) 15.18,38.43 

-0.59X10 

2.863 


( . 1/ — / > >22, 4 1.83. 71) 

-0.59527 

2.803 


(. 1 / - nr 4lA2A7 

-0.59414 

2.864 


( .4/ I * ),3$,42.43 

-0.5904X 

2. 8(54 


( .4/ — T I ) 38.39, 41 .42 

-0.58915 

2.864 


(Al — / / ) 18 , 19,42.43 

-0.5X037 

2.864 


M/ - r* 

-0.5X030 

2.864 


( A i — 1 i > 52-08 + ( -4/ — A i )o . r >-88 

-0.351X2 

2.XX2 

6.25 

(.4/ - rors 

-0.59975 

2.865 

6.94 

(.4/ — Ti) 22,4 1.64 .(>(>,72 

-0.59194 

rag 


( 4/ — I ♦ )ir>J 8.38,43. 4«'i 

-0.59098 

1 

7.11 

(*4/ I *) 14.19.39,42 

-0.593X3 

2.868 


( -4/ — 1 1 ) 14.42, 43.07 

-0.5935X 

2.868 


( .4/ / 0 14. 4 1.43.70 

-0.5X945 

2.868 


( *4/ 7 / )* 5.19.38,42 

-0.5X502 

2.869 


( -4/ ^ / )37,42. I3.W 

-0.5X500 

2.869 


( -4/ 1 ' >17,42.85.07 

-0.5X52X 

2.869 


( 4/ 7 i ).38.4 1 .13.4(1 

-0.5X400 

2.869 


( *4/ 7 / 1.38.39, 42. 13 

-0.57X47 

2.869 


( 4/ 7 t >18.38,42.48 

-0.57X03 

2.869 
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( 'on figuration 
( - 1/ — 77 )ir>. 18.20,38.40.43 

(-1/ — I i >38. 10.11.13,40.48 
(Al — / i )20-22, 40.02. 08 
( Al — Ti )i 7. 19.4 1, 43.03,07 

( .4/ — 7 / )i7,-n-n.<i r > 

( .1/ 7 / ) 14,19.39.12.09 

( ^ ^ * )|1. 19.22, 39, 12 

( 1 * ):^7. 12,43.02,70 

( Al 1 i ) 1. 1,22,37.39.12 

(Al 1 i )]. 1.37,39. 4 2. 40 

( Al / 038.39.12,43 

(Al ^ O^IT. 19, 12.05,07 

( Al 7 t ).*9,42.43,45.4o 

(Al 7 0:^,41.42.13.40 

(Al — 7 18.38,41, 43, 40 

(•4/ — 7 /) 17/23, 39.17, 01.00.72 

( Al — / > )i3.13. 18.20.38.40.41. 13 

( Al 7 / );^ 4 j I .00.1 >9 
(.4/ 7 / )l3. 18,38, 41 .13.00 

(.4/ — 7 / )42.44.4M8,05.07m71 

( A l ^ ^13.21 .38.10.03.71 

(-4/ 7 / >lV. 19.22,37.39,07 

(i4/ 75 )i3,l 7.19, 12.45,00 

7'/) 14 .18.37,39,47,00 
(.4/ — Ti )l9,20,37,39.4M4.0U,U7 
( *4 / — T / ) 15.1 9,38 .4 2 ,U 1 ,( >4 .( 49.72 

(Al Ti )|j8.41 — 43.4U 

(*4/ — Ti )imj 7.21, 03.07, 71 

( *4/ 7 / ).37 , 4 2 ,43 ,4 5 ,47 

(,4/ — I i )n. 17,18,21.42,00 
I (Al 7 037,39,4 1—43.4U 

(.4/ — T i ) 18,19.22,23, 42. 43, 10. 17 
(-4/ — 7 041-44,(i- r >-(i8 

( .4/ — T l )i3 42 

(.4/ — T 0 13,1 5.18,20,38.40, 41 ,43,10 
(.4/ — Ti ) 13,19.22-24,38.43.02.67 

(Ai-Tiy^ 

(Al — 77 )a8— 40,42 — 14,40 — 18 

(.4/ — 7 i ) 14,1 8,22,38,42,40.02.00.70 


A// (oV/alom’ 
-0.59 151 
-OASIS-i 
-0.5X245 
-0.57174 
-0.57094 
-0.5XX92 
-0.5X105 
-0.5S05S 
-0.5X02X 
-0.5799X 
-0.57470 
-0.57412 
-0.57204 
-0.50017 
-0.5X444 
-0.5X159 
0.5 XI OS 
-0.57594 
-0.57552 
-0.57402 
-0.57245 
-0.57197 
-0.57157 
-0.50X02 
-0.50079 
-0.50002 
-0.50170 
-0.50110 
-0.50050 
-0.55X17 
-0.55792 
-0.55445 
-0.55220 
-0.5X400 
-0.5X000 
-0.577X9 
-0.50041 
-0.55910 
-0.541X4 
-0.54444 


NASA TM-1 13121 



• r Ti 

( 'on figuration 

A// (oV/atoin) 

« (A) 

12.96 

(.1/ I 

-0.570S7 

2.8X5 


(.1/ 1 0i:},ir,.i8.38.-i:u«».(iJ> 

-0.5 70S I 

2.885 


(.1/ I i )|8,:i8.-u.*ni.ni.<ici,<in 

-0.56697 

2.884 



-0.56550 

2. SSI 


(.1/ 7 '0 19 . 21 .38. 11.12. K). 17 

-0.50090 

2.8X4 


(.1/ 7 / ) 1 8.37.39,4 2 .4 5 .4 7.< ill 

-0.55261 

2.885 


(.4/ 1 t )ji 18.22.37,38.10, 17 

-0.55145 

2.SS5 


(A/ 7^ ) 18.38,4 1-43,46 

-0.54795 

2.885 


^ ^ * ) 38 ,39, 41 — 43.4 5 .46 

-0.5 1300 

2.XX0 


(-4/ ^').38. 39,41 -43.46,47 

-0.51099 

2.XX0 

14.89 

( 4/ — 1 i ) | 3 . 1 5,18.20.2 1 ,38.40.4 1 .43.40 

-0.57145 

2. 885 


(-4/ — 7 /)37.30,42,44.45,47.02.04,t>- r >,07, 70.72 

-0.55922 

2. 895 


(.1/ — 7 /)ir,. 18.10.38, 41,42.44.47.<v3.(i9 

-0.55710 

2.887 


( -4 / — / / ) 1 5. 1 9,23,37,39,4 1 .43.4 5.47,6 1 .05 ,09 

-0.52729 

2.S96 


( .4/ — 1 l )37-48 + ( A' *49-60 — -4/ ( jl_72 ) 

-0.52156 

2.905 

1 LSI 

(.4/ — 7 i)] { 15.17,22.37, 43,02. 71 

-0.56500 

2. XXX 


( -4/ I i 4 1 5.1 8.23,38.43,01 ,03.00 

-0.55915 

2.SSS 


(Al— 1 i) 14 17,1 9,22.37.39.4 2 .4 5 

-0.55SS7 

2.SSS 


(-4/ l * ) 1 :5. i 7.18.22.38.41 .42.45 

-0.547X1 

2.889 


( *4 / 7 0l8.iO.22, 23.41.42, 45.40 

-0.54140 

2.890 


( -4 / — 7/')| 3,m. it, 18,37.38.41. 42 

-0.556S9 

2.890 


(.1/ — / /) 13.14.18,22,37.38,42.40 

-0.55557 

2.89 1 

16.67 

( J / ? / ) 1 5 1 8.23,38, 4 1 .13 .01 .00. 7 1 

- 0,55407 

2.895 


(.4/ 1 Oh. 19,21, 37,39.41. 42.40. 17 

-0.54S40 

2.895 


( .4/ '7 0.38.41,43,46,6! ,63,66.69.71 

-0.54782 

2.895 


( * T * ) H. 1 7,22 .39,4 2 ,4 7 .6 1 .05 .09 

-0.54745 

2.895 


(*4/ — I i ) 13 .23.38,39,42,45.40.01.71 

-0.54405 

2.894 


(.4/ 7’/)| f>ly .23,37.41,45.02,00.70 

-0.54747 

2.895 


7 / ) ;w ii ,42, 15,46.65,66,69.70 

-0.54216 

2. 895 


7’/)i 8 3 7 _39 42 ,45~47,W 

-0.52517 

2.896 


( *4/ 7 Ooi— 03,65-07.09-71 

-0.5 12 IS 

2.897 

IS. 52 

7 0 13. 18.21,38,39.41, (J2.07.09.70 

-0.54459 

2.898 


7 Oh ,17 .jg ,42.01.03,05,07,09.71 

-0.55964 

2. 899 


7 0 14,17. !9,22. 4 2. 01, 03,00 .60. 71 

-0.5595S 

2.899 


(-4/ 7 > ) 13,17,38,41, -43,45,47.0(1.70 

-0.52875 

2.900 


(.4/ 7 /) 18 , 19 . 38 , 39 , 42 . 13 , 03 . 05 .(^ 1.71 

-0. 5258 7 

2.900 

IS. 75 

7 / ) 13 . 18,38 

-0.50750 

2.X90 


( 7 ; ) 13.17,38 

-0.54137 

2.X 90 


* r? 

* — 

£ 

1 
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tc 


Figure 1 . — Schematic illustration of the Heusler phase and Its relation to B2 
compounds. 



X/Ni 

Figure 2. — Lattice parameter of non-stoichiometric NiAl alloys as a function of 
Ni concentration, normalized to their stoichiometric values. The solid 
squares denote results from different investigators (see ref. 1). The lines 
denotes the BFS predictions. 
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(a) 



Figure 3. — (a) The computational cell used to generate the catalogue shown 
in appendix B. (b - c) Two samples of the configurations of NiAI-Ti alloys used 
in this calculation, corresponding to Xjj = 1 1 .1 1 and = 1 6.67. 
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Figure 4. — (a) Energy of formation (in eV/atom) of the cells listed in 
Appendix 2 for 0< x-pj < 50. (b) A subset of the configurations listed 
in Appendix 2, for the range 0 < x-n < 25. The circles denote those 
configurations characterized by Heusler ordering, while the solid 
squares include a variety of short-range order patterns as well as 
disordered states. 
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Figure 5. — Configurations used for the calculation of Ti site preference 
in NiAl alloys, (a) A Ti atom in an Al site, (b) a Ti atom in a Ni site, 
with the displaced Ni atom occupying a site in the Al sublattice at 
nearest-neighbor distance from the Ti atom, (c) same as before but 
with the Ni atom located at a different site in the lattice, at a distance 
greater than next-nearest-neighbor distance. In all cases, A (B) 
indicates an atom A in a B site, while A(B)C indicates an atom A in a 
B site with the displaced B atom in a C site, the B and C sublattices 
correspond to that of an ordered B2 structure. 
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Figure 6. — Energies of formation (in eV/atom) for 72 atom cells containing 36 
Ni atoms, 36-x Al atoms and x Ti atoms (x = 1 , 2, 5). The different energy 
states correspond to different substitutional defect schemes, as indicated 
in Fig. 5. 


• A 



Figure 7. — Schematic representation of Heusler and Kneen ternary phases, (a) Heusler phase A 2 BC. (b) Kneen 
phase A 2 BC. 
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(b) 


Figure 9. — Final structures of a Monte Cailo/Metropolis/BFS 
simulation on a 1 024 atom cell of a Ni-25AI-25Ti alloy. Both, 
initially random, states are obtained by lowering the temp- 
erature in different ways: Fig. 9a shows the final geometry 
for a rapid cooling process, where the final temperature is 
reached by 'freezing' the initial, high temperature, state. 

Fig. 9.b is obtained by slowly lowering the temperature in 
equal 1 00 K, temperature intervals until the final temperature 
is reached ('cascade' process). The inset includes the 
coordination matrix for the final state of the simulation (see 
text): the element ij in this matrix indicates the probability 
that an atom / has an atom of species j as a nearest- 
neighbor. The labels /, j take the values 1 , 2 and 3, corre- 
sponding to Ni, Al and Ti. 
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Figure 1 0. — Final geometry for a cascade simulation (see text) on a 
Ni-49AI-1Ti cell of 1024 atoms. An expanded view of the cell is 
also shown. 
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Figure 11 . — Results of three consecutive temperature cascades on a Ni-45AI-5Ti 
1024-atom computational cell. The fourth column corresponds to a cascade 
calculation on a larger (4608 atoms) cell. 
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Figure 1 2.— Final geometry for a double cascade process on a Ni-40AI-10Ti 
1 024-atom cell. The first expanded cell shows the results of the first cascade 
process, where the cell is slowly cooled at equal temperature intervals from 
an arbitrary high temperature. Then this cell is re-heated and slowly cooled 
again, reaching the final state shown in the second cell. 
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Figure 1 3.— (a) Bright-field TEM image of the micro- 
structure of the NM7A1-3T1 alloy and (b) corres- 
ponding <1 1 0> SADP. The only features observed 
in the alloy were an occasional <1 00> dislocation 
marked "D" on (a). 
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Figure 14. — (a) Bright-field and (b) dark-field TEM image of the Ni-45AI-5Ti alloy showing precipitation of fine 
Heusler particles and (c) corresponding SADP and (d) an indexed, simulated pattern. 
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Figure 1 5. — Dark-field TEM image of the Ni-43AI-7Ti alloy showing 
dense precipitation of rectangular-shaped Heusler precipitates. 
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Figure 1 6. — Summary of final states for cascade processes on 
Ni-(50-x)AI-xTi alloys (x = 1 , 5, 6, 7, 8, 1 0, 25) for 1024 atoms cells. 
The Heusler alloy Ni-25AI-25Ti results correspond to a first cascade 
and a third cascade process, highlighting the stability of the heusler 
phase (see text). 
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